Exceptionally large room-temperature ferroelectric polarization in the novel PbNiO,3 

multiferroic oxide. 
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We present a study based on several advanced First-Principles methods, of the recently synthe- 
sized PbNiC>3 [J. Am. Chem. Soc 133, 16920 (2011)], a rhombohedral antiferromagnetic insulator 
which crystallizes in the highly distorted R3c crystal structure. We find this compound electrically 
polarized, with a very large electric polarization of ~ 100 /xC/cm 2 , thus even exceeding the polar- 
ization of well-known BiFeC>3. PbNiC>3 is a proper ferroelectric, with polarization driven by large 
Pb-O polar displacements along the [111] direction. Contrarily to naive expectations, a definite 
ionic charge of 4+ for Pb ion can not be assigned, and in fact the large Pb 6s-0 2p hybridization 
drives the ferroelectric distortion through a lone-pair mechanism similar to that of other Pb- and 
Bi-based multiferroics. 

PACS numbers: 77.80.-e, 77.84.-s, 77.22. Ej, 71.20.-b 
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Introduction Multiferroics are materials in which differ- 
ent ferroic orders such as ferromagnetism, ferroelectricity 
and/or ferroelasticity may coexist in one single phase 
In the last few years, there has been a tremendous boom 
of interest in these materials, due to the potential appli- 
cations in memory devices or in novel type of magnetic 
switching in magnctoclcctric multiferroics, based on the 
cross-coupling between ferroelectric and magnetic chan- 
nels Furthermore, these materials offer a rich and 
fascinating playground for the complex physical mecha- 
nisms underlying the processes involved in the observed 
properties. [B} It is obvious that the search and the pre- 
diction by material design of new multiferroics is of great 
importance for both fundamental physics and technolog- 
ical applications 0, 0] • 

Relatively few multiferroics have been identified so 
far. 0] Without doubt, the most studied and well char- 
acterized multiferroic is BiFeC>3 [9]. It crystallizes in the 
polar space group R3c (No. 161, point group Cs v ), and 
is predicted to have a G-type antiferromagnetic (AFM) 
alignment of the Fc spins @, Ulj . The R3c symme- 
try corresponds to the so-called LiNb03-type structure, 
which can be viewed as a highly distorted double per- 
ovskite with rhombohedral symmetry. The primitive 
unit cell contains two formula units (10-atoms), arising 
from counterrotations of neighboring O octahedra about 
the [111] axis (see FigfJ}. The crystal symmetry al- 
lows the presence of a spontanous polarization along the 
[111] direction, which arise from the relative displace- 
ment of the Bi sublattices with respect to the FeOg octa- 
hedra cages along [111]. The origin of the large sponta- 
neous polarization of Bi +3 Fe +3 0~ 2 3 , ~ 90 /LtC/cm 2 0, 
has been explained by first principles calculations within 
the density functional theory plus Hubbard-U approach 



(DFT+U 

(mtp) Ha 




and the "modern theory of polarization" 
, in terms of the lone-pair electrons present 



at the Bi sites which are ultimately responsible for the 



large displacement along the [111] direction Q. 

Inspired by the recent reports of Inaguma et al. 
on the synthesis of a new antiferromagnetically or- 
dered compound with LiNb03-type structure, such as 
PbNiOa [12], [l3| , we explore here the possibility of mul- 
tiferroic behavior in PbNi03. We first summarize the 
experimental findings. Inaguma and coworkers synthe- 
sized two high-pressure polymorphs of PbNi03 with a 
(A) perovskite-type structure and (B) the LiNbC>3-type 
structure [Hill; the latter (hereafter called L-PbNi0 3 ) 
is thermodynamically more stable than the perovskite- 
type one at ambient pressure. With respect to the 
orthorhombic structure, in the acentric rhombohedral 
LiNb0 3 -type (R3c) structure (Fig. [TJ Pb and O atoms 
are displaced against each other along the threefold [111] 
axis leading to a large distortion of the PbOg and NiOe 
octahedra (see Fig[TJ) . The Pb atom is coordinated by 
six oxygens at 2.10 and 2.25 A, while the Ni-0 bond 
distance splits into two subgroups (2.07 and 2.11 A). 
Magnetic susceptibility and resistivity measurements in- 
dicate that L-PbNi03 undergoes an AFM transition at 
Tn=205°K, and exhibits semiconducting behavior. The 
AFM ordering in the acentric crystal structure suggests 
possible multiferroic behavior. 

In the following we show that L-PbNi03 is a new room- 
temperature multiferroic, with an exceptionally large po- 
larization (P) of about 100 /iC/cm 2 , which is the highest 
polarization ever predicted for any bulk material so far, 
and ?» 10% larger than that of BiFe0 3 @. 

Computational details We performed DFT based cal- 
culations using the Vienna Ah initio Simulation Pack- 
age (VASP) 11611 based on the projector-augmented-wave 
method [13, EH, withing the Perdew, Burke and Ernz- 
erhof (PBE) parameterization scheme (l5| for the gener- 
alized gradient approximation (GGA). To overcome the 
deficiencies of standard exchange-correlation approxima- 
tions for localized Ni d states, we made use of three 
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FIG. 1: (Color online) Schematic view of the LiNb03-type 
PbNiC>3 within space group R3c structure built up from two 
cubic perovskite unit cells, the black (large), gray (medium) 
and red (small) spheres denotes Pb, Ni and O atoms, respec- 
tively. With respect to the centrosymmetric (R3c) paraelec- 
tric phase (see text) the cations are displaced along the [111] 
axis relative to the anions (the arrows represent the displace- 
ment vectors), and the oxygen octahedra slightly rotate with 
alternating sense around the [111] direction [Q. The crystal 
structure was drawn using the program VESTA[l4l]. 



beyond-DFT approaches: (i) Dudarev's GGA+U 0, 
|20| . using U=4.6 eV in accordance with constrained DFT 
calculation of Refill! (ii) The renowned Heyd-Scuseria- 
Ernzerhof (HSE) screened hybrid functional [ID], which 
has been shown to give an excellent account of materi- 
als properties for magnetic multiferroics 23] . (iii) The 
recently introduced variational pseudo self-interaction- 
correction method VPSICJIJ] (adopting the HSE opti- 
mized structure), implemented within plane- wave basis 
set and ultrasoft pseudopotential scheme in the PWSIC 
code. The cutoff energy was set as 600 eV and a 8x8x8 
Monkhorst-Pack grid of k points was used. In HSE, the 
fraction (1/4) of Fock exchange was sampled using the 
twofold reduced k-point grid, to reduce the computa- 
tional load. The lattice parameters and atomic positions 
were relaxed (at GGA+U and HSE level) until the total 
energy changed by less than 10~ 5 eV per unit cell and 
the residual force was smaller than 0.01 eV/A. 

Magnetic, structural and electronic properties The G- 
type antiferromagnetic configuration (AFM-G), in which 
each Ni ion is surrounded by six nearest-neighbors Ni 
with opposite spin direction, is the most favourable 
spin configuration with respect to ferromagnetic (by 50 
meV/f.u.), AFM type-A (by 33 meV/f.u.) and AFM 
type-C (by 17 meV/f.u.) orderings. 

The resulting calculated magnetic Ni magnetic mo- 
ment (mm) is about 1.7 hb at both GGA+U and HSE 
level. The optimized GGA, GGA+U and HSE structural 



TABLE I: Structural parameters, magnetic moment and elec- 
tronic band gap for AFM-G configuration of L-PbNi03 within 
GGA, GGA+U and HSE method, compared with the avail- 
able experimental data, a, c are the lattice parameters in 
the hexagonal setting, whereas x, y, and z are the internal 
atomic positions of Pb and O. Ni ions sit in the (6a) (0,0,0) 
positions. Bond lengths of Pb-0 and Ni-O are also reported. 
Ni — O — Ni (°) is the Ni-O-Ni angle, m indicates the mag- 
netic moment of Ni. 





GGA 


GGA+U 


HSE 


Expt. 


a (A) 


5.442 


5.430 


5.359 


5.363 


c(A) 


13.737 


14.353 


14.209 


14.090 




0.2226 


0.2885 


0.2898 


0.2864 


xo 


0.0605 


0.0532 


0.0468 


0.0487 


yo 


0.2598 


0.3543 


0.3579 


0.3657 


zo 


0.0931 


0.0702 


0.0687 


0.0668 


Ni-0 (A) 


2.051 


2.128 


2.119 


2.109 




1.905 


2.060 


2.053 


2.071 


Pb-0 (A) 


2.395 


2.310 


2.258 


2.246 




2.238 


2.174 


2.121 


2.104 


Ni^~0^ Ni (°) 


147.8 


140.6 


138.1 


136.7 


mm (a»b) 


0.97 


1.67 


1.69 




gap (eV) 


Metallic 


0.37 


1.18 





parameters for the AFM-G phase are listed in TableQ] As 
expected, GGA does not reproduce well the experimental 
values, whereas HSE and, to a lesser extent GGA+U, de- 
liver numbers in good agreement with experiment, with 
relative errors of 1.5-2 % (GGA+U) and < 1.0% (HSE). 

The drawbacks of conventional GGA are more dra- 
matic for the electronic properties. The GGA scheme 
finds a metallic character, in clear disagreement with ex- 
periment (l2| . Conversely, the inclusion of either the on- 
site U or the Fock exchange cures this limitation and 
correctly predicts an insulating state with an energy gap 
of - 0.4 eV (GGA+U) and - 1.2 eV (HSE), as inferred 
from the density of states (DOS) shown in Figj2](the VP- 
SIC gap is 0.9 eV, see Supplementary Materials (SM)). 
The gap opens between occupied Ni d and O p states, and 
the lowest unoccupied electronic states that are made up 
by a mixture of Pb s and O s states. 

Spin polarization does not affect Pb and O DOS, 
whereas the Ni atoms display a magnetic moment of 
1.69 /ib- In hexagonal coordinates (with z parallel to the 
[111] direction) the Ni d states are distributed into two 
doublets, (d xy , d x 2_ y 2) and (d xz , d yz ), which contribute 
to the net magnetization by 0.65 and 1.65 /ib, respec- 
tively, and a spin-degenerate d z 2 singlet (see SM for an 
Zm-projected DOS). The broad Ni 3d valence manifold is 
strongly hybridized with the O 2p states, and lies in the 
energy range between -6 eV and the Fermi level. 

The most important feature of the DOS is at the bot- 
tom of the valence region, between -9 eV and -6 eV: a 
large region dominated by highly- hybridized Pb 6 s and 
O 2p states. The presence of a substantial amount of 
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charge (about 1 electron, see SM) within the Pb 6s or- 
bital is at odds with the nominal Pb 4+ valence originally 
assumed in the interpretation of X-ray photoemission 
spectroscopy (XPS) results 12|, [l3j . In realistic calcula- 
tions, substantial deviation between the static charges 
and the nonimal ionic charges (due to large hybridiza- 
tion effects) is the rule rather than the exception, thus it 
certainly does not come as a surprise. In this case, how- 
ever, this deviation is worthy to be emphasized, since, 
as we will show below, it is strictly connected to the fer- 
roelectric instability. To highlight the role of Pb Qs-0 
2p hybridization, in Figj2]we show the DOS for the cen- 
trosymmetric, non-polar reference structure R3c as well 
as the polar structure: the major difference is indeed the 
2 eV downshift and the broadening (more than a factor- 
2 bandwith expansion) of Pb 6s-0 2p spectral weight 
occuring along with the centrosymmetric-to-ferroelectric 
transformation. This so-called stereochemical activity of 
the A-site cation [3, [29| is quite consistent with what 
occurs in Bi-based (BiFe0 3 [l, BiMnOalU, 0]) and Pb- 
based (PbTiO3[30]) compounds, usually labelled as lone- 
pair ferroelectrics. (see SM for a more detailed descrip- 
tion of the bonding picture in PbNiOs and a one-to-one 
comparison with the BiFeC>3 case). This is also reflected 
in the structural distortions which predominantly involve 
counterdisplacements of Pb and O atoms along the [111] 
direction (see FigHJ). 

In the following we explore the ferroelectric properties 
of L-PbNiC>3 by evaluating the spontaneous macroscopic 
polarization P. 

Analysis of ferroelectric polarization 

By pseudosymmetry analysis [31] we determine a par- 
ent centrosymmetric structure (R3c) with minimal su- 
pergroup symmetries, from which the observed non- 
centrosymmetric structure can be reached through a con- 
tinuous structural distortion. By definition, the polariza- 
tion is zero (modulo a polarization quantum 27. 32|) in 
this paraelectric (PE) phase (A = 0), thus it represents 
a viable reference state for the evaluation of the polar- 
ization in the polar structure (A = 1). Here A = 1 is 
the amplitude of the polar distortion which progressively 
transforms the R3c phase into the R3c one. For consis- 
tency, we assume for the R3c phase the same volume and 
rhombohedral angle of the polar structure. 

In the PE phase, Pb and O atoms lie in (111) PbC>3 
planes. Symmetry mode analysis shows that the A = 
[0,1] transformation (showed by arrows in Fig.l) is al- 
most exclusively composed by in-phase [lll]-parallel rel- 
ative shift of Pb and O atoms (i.e. a nearly-exclusive 
r2_-mode transformation), with some very residual O- 
rotation contribution. The Pb atoms are displaced from 
their positions by 0.56 A and the O atoms by 0.23 A, 
while Ni atoms are almost unchanged. This transforma- 
tion produces an alternate shrinking and dilatation of 
Pb-0 distances along [111] of about 0.8 A, whereas Ni 
(111) planes are left off-centered by 0.17 A with respect 
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FIG. 2: (Color online) HSE calculated total and site- 
decomposed density of states for the AFM-G type PbNiOs in 
the acentric R3c (full lines) and centric R3c (shadow) phases. 
The vertical line denotes to the Fermi level. 



to their adjacent PbC>3 planes. The total polarization 
P can be split as the sum of ionic Pi on and electronic 
contribution P e ; e : the former is the dipole of the ion- 
core charges; the latter is obtained by Berry phase ap- 
proach within the MTP 



27, 321. HSE and VPSIC de- 



scribe both polar and nonpolar structures as insulating, 
whereas within GGA+U a large U=7.6 eV is required to 
open a gap in the PE phase. The HSE-calculatcd DOS 
shown in Figj2]indicates that the value of the gap is essen- 
tially the same for both acentric R3c and nonpolar R3c 
phases. Fig. [3] displays total energy and polarization P to t 
as a function of the polar distortion A at both HSE and 
PBE+U level. Both methods deliver the same outcome. 
As expected the polar structure is more stable than the 
non-polar one by 0.6 eV/f.u (0.7 eV/f.u. according to 
VPSIC). Both Pj on and P c i grow monotonically with 
the polar distortion, giving rise to a total polarization 
Ptot of - 100 fiC/cm 2 (98.5 fiC/cm 2 according to VP- 
SIC), thus 10% larger than that of BiFe0 3 [a El and 
significantly enhanched with respect to that of LiNbOs 
(80 ^C/cm 2 )[IH and ZnSn0 3 (57 fiC/cm 2 )^. 

We have also evaluated the polarization using the 
approximated expression Ptot = XL Ai?^ Z*, where 
ARi are atomic displacements (from centrosymmetric 
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FIG. 3: (Color online) Total energy profile -E, total polariza- 
tion Ptot as a function of the polar distortion A, from the 
A=0 PE phase to the A=l ferroelectric (FE) one. Both HSE 
(dashed lines, dark) and GGA+U (dotted lines, gray) data 
are displayed. 



TABLE II: Born effective charges for Pb and Ni in PbNi0 3 
(with ferroelectric and paraelectric phases) computed via 
linear-response formalism within the GGA+U method. 





Pb 


Ni 


O 


R3c (FE) 


4.40 


1.97 


-2.12 


R3c (PE) 


4.41 


1.97 


-2.12 



and comments. A.S. thanks J.M. Perez-Mato for use- 
ful discussions. The calculations have been performed 
on the Vienna Scientific Cluster (VSC) and, partially, in 
CASPUR Supercomputing Center in Rome. 



to polar structure structure) and Z* are Born effective 
charges (BEC) calculated via linear-response formalism 
in GGA+U. We obtain a value (93.8 ^C/cm 2 ) substan- 
tially similar to the evaluation based on the exact Berry 
phase formula. We found a nearly isotropic BEC tensor 
and dynamical charges similar in both PE and FE phases 
(see TabHJ). For Pb we obtain an anomalous BEC value 
of 4.4 (as compared to the static charge, « 2.15), identical 
to the Bi BEC calculated in BiFeOs, coherently with the 
role of this ion in guiding the ferroelectric transition^. 

Summary Using several beyond-LDA methods 
(LDA+U, HSE, and VPSIC) and the MTP approach for 
the determination of polarization properties, we report 
the presence of a large spontaneous electric polarization 
~ 100 /iC/cm 2 in rhombohedral PbNiOa associated to 
the structural transformation from the centrosymmetric 
R3c to the polar R3c symmetry. The microscopic 
analysis shows that this proper ferroelectric instability 
is dominated by huge in-phase (i.e. T-mode) relative 
displacements of Pb and on-top O atoms along the [111] 
direction, in turn associated to large Pb s-0 p band 
rehybridization at the bottom of valence band manifold. 
This ferroelectric mechanism can be then assimilated 
to other Bi-based and Pb-based ferroelectrics, whereas 
the naive interpretation of Pb as an inactive +4-charged 
ion would be totally misleading. This material may 
be prototype of a new class of Ni-based rhombohe- 
dral multiferroics, which take advantage of the stable 
rhombohedral symmetry to develop large ferroelectric 
displacements along the [111] axis, and still mantain 
strong magnetic coupling in the Ni sublattice. 
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